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Abstract: The three 3-brane system with both positive or negative tension is stud- 
ied in a low energy regime by using gradient expansion method. The effective equa- 
tions of motion on the brane is derived and in particular we examine, in the first 
order, the radion effective lagrangian for this system. In this case, we show the so- 
lution of the modified Friedmann equation with dark radiation on the middle brane 
and the other 3-branes by direct elimination of the radion fields and Weyl scaling of 
the metric on the branes. We also derived the scalar-tensor gravity on the branes. 
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1. Introduction 

One of the long standing problem in particle physics and gravitational theories is how 
to understand quantum theory of gravity. Nowadays, the only possible candidate for 
this theory is the superstring theory Jl|. Interestingly, this theory predicts the exis- 
tence of the extra dimensions. In order to reconcile this prediction with our observed 
four dimensional universe, we need a mechanism to compactify the extra dimensions. 
In particular the setup of heterotic M-theory and its compactification down to five 
dimensions || leads to a well motivated five dimensional brane world scenario, which 
can be used to study its consequences in particle physics and cosmology. 

Randall and Sundrum (RS) proposed two similar but distinct phenomenological 
brane world scenarios ||, |J. The first scenario is composed of two branes of the op- 
posite tension, namely RS I. This scenario is the five dimensional space-time which 
all matter fields are assumed to be confined on branes at fixed points of the S 1 /Z2 
orbifold so that the bulk is described by pure Einstein gravity with a negative cos- 
mo logical constant. On the other hand, the second scenario has a single brane with a 
positive tension, namely RS II. The fifth dimension is infinite but still Z 2 symmetry 
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is imposed. In both scenarios the existence of the branes and the bulk cosmological 
constant makes the bulk geometry curved, or warped. 

Furthermore, the brane world models are expected to shed some light onto not 
only quantum gravity and unification issues, but also cosmological issues such as the 
cosmological constant or dark energy problem, or even particle physics ones such as 
the hierarchy problem. There is also a brane world alternative to the standard big 
bang plus inflation scenario 0, E|. This scenario consist of a five dimensional bulk 
bounded by two branes, which are as usual located at the fixed point of the S 1 /Z2 
orbifold. Although many brane world models might be over simplified, they should 
help in learning and understanding the properties of an effective theory derived from 
some action in space-time dimensions d>4 0. 

There are two approaches to obtain the effective Einstein equations on the brane 
in the context of Randall- Sundrum scenarios, namely covariant curvature formula- 
tion and gradient expansion method. In the covariant formulation ||, the effective 
Einstein equation can be obtained by projecting the covariant five dimensional Ein- 
stein equations on to the brane. The resulting projected equations are modified with 
respect to general relativity due to the presence of a local quadratic term in the 
sources and to the presence of a non-local term which is the projection of the five di- 
mensional Weyl tensors. This last term carries information of the bulk gravitational 
field on the brane and its contribution is of fundamental importance as it might be 
relevant even at low energy [H| [11], [12 . 



The main difficulty in understanding the contribution of the projected Weyl 
tensor to the effective theory on the brane is in its non-local character. The equation 
for the projected Weyl tensor on the brane are not closed so that solving the full five 
dimensional equation of motion is necessary. 

On the other hand, gradient expansion method gives a way out of this problem. 
This method first proposed by Kanno and Soda [ll|. The main idea is to treat the 
issue perturbatively, defining a low energy regime in which the energy density on 
the brane is kept small with respect to bulk vacuum energy density. The perturba- 
tion parameter is defined as the ratio between these two energy densities. The five 
dimensional equations of motion can be solved at different orders in the perturba- 
tion parameter. This method allows in principle to derive the effective equations of 
motion on the brane at each order. 

In this paper, we generalize the results obtained for a two 3-brane scenario (RS 
I) to a multibrane scenario [|l4], 15, 16 1 with the set up of the following. We consider 
the A, B, C 3-branes system, with A and B branes are placed at the fixed points 
of the orbifold whereas C brane is put between the two fixed points. We derive the 
effective equations of motion in this scenario. By performing a perturbative expan- 
sion of the metric, an expansion of the extrinsic curvature tensor and Weyl tensor 
are considered. The four basic equations for the five dimensional evolution equations 
and junction conditions are then solved at different order in the expansion parameter. 
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The parameter of expansion is determined as in the Anti-deSitter (AdS) scenario. 
There is a constant scale, namely the AdS curvature scale, to which quantities can 
be compared. 

The paper is organized as follows. In section 2, we explain the general set up for 
the multibrane system. The junction conditions for this system is derived for each 
brane by using geometrical approach. In section 3, we give the basic formulation 
of the gradient expansion method by solving five dimensional equations of motion 
and imposing the Dirichlet boundary condition at the brane position. We give the 
basic formulations of the extrinsic curvature at the zeroth order and the first order 
expansions. The first order solution of the effective equations of motion on the brane 
is derived in section 4. In section 5, we address our result into the scalar-tensor 
theory. We conclude our results in section 6. 



2. The general setup and background solution 

In this section we discuss three 3-brane embedded in a five-dimensional space-time. 
This model is a straightforward extension of the RS I model. In more detail, we have 
three parallel 3-branes in an AdS 5 space with negative cosmo logical constant. The 
fifth dimension has the geometry of an orbifold S 1 /Z 2 and the branes are located at 
y = (A-brane with tension a a), y = yA (C-brane with tension ac) and y = ys (B- 
brane with tension <tb). Note that A-brane and B-brane are placed at the fixed points 
of the orbifold S l / Z 2 but C-brane is not. The region between two brane, namely 
region A [0, jja] and region B [yA, Vb], is characterized by two different values in two 
different slice of AdS$ curvature scales i.e., Ia and Ib, respectively. Furthermore the 
metrics on the three branes are all connected by an conformal transformation, so in 
principle it is enough to derive the four- dimensional effective equations of motion on 
one of the branes. In the following we derive the four-dimensional Einstein equations 
of motion on each branes. The action that describes this configuration is 

S = ^ J <?X^-g (n+^j+ J dW-g'-^i-Vi + ^matter), (2-1) 

where TZ, k and g are the scalar curvature, the gravitational constant in 5-dimensions, 
respectively. The metric g l - brane is the induced metric on the branes and Oi their 
tensions. Notice that we have allowed the curvature scales to have different values 
to either side of the middle brane, represented by I a and Ib, which leads to different 
metrics in the two bulk regions. The coordinate system chosen is 

ds 2 = eWv^dy 2 + g^(y, x)dx t "dx v . (2.2) 
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The proper distance between two branes with fixed x coordinates can be written 



as 



d A (x) = / dy'e^ & , (2.3) 
Jo 
rVB 

d B (x) = / dy'e^'^ . (2.4) 
JyA 

And the total proper distance is given by 

rVB 

d t ot(x) = d A (x) + d B (x) = / dtfeMrt . (2.5) 

Jo 

The Einstein equations that arise from this coordinate system are: 

e" W),y - {e-+K){e-*K») + < 4 >i£ - D,D»cf> - D^W* 
= -j^S; + (±5;a A + If - l -5^ e-%y) 

+y (fac + n v - \ s ; fc ) - va) 

+k 2 (^S;a B + 7f - i^f B ) e-^(y - y B ) , (2.6) 



and 



e -*( e -<l>K) ty - (e^K^e^K^) - D a D a (f) - D a <j>D a <j> 

4 ^ (-4a A + T A ) e-*8{y) - y (-4<r c + f c ) e~H{y - y A ) 



I 2 

,2 



-y (-4a B + f B ) e-^(y - y B ) , (2.7) 
and 

A,(e-^)-Zy e -^)=0, (2.8) 

where the appropriate AcLSs curvature scale is chosen for each region. is the 

curvature on the brane and denotes the covariant derivative with respect to the 
metric g^ u . When the metric changes between regions of the bulk, it is necessary 
to include a term in the boundary action that depends on the trace of the extrinsic 
curvature. 

The extrinsic curvature tensor is defined by K^ v = —g^ y /2. For all three branes, 
the junction conditions read 

e" [K; - S;K} \ y=0 = ^ (-a A S; + Tf) , (2.9) 

e" K - <M - \y=y- = ^ ("^ + if) , (2.10) 

e" K - I™ = " Y ("^ + ^ ) > ( 2 - n ) 
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where we have used that A-brane and B-brane follow a Z2 symmetry but C-brane 
is no Z2 symmetry. Here, we assume that the direction of the normal vector field to 
a brane is chosen to be the same all way through the bulk for all the three brane. 
Decompose the extrinsic curvature into the traceless part and the trace part 



e~*K, u = + - A g„ v Q , Q = -e~*^- log 



then, we obtain the basic equations; 



-<f>spu 



QK 
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DAD" 



+- A K (D a D c 



D a 4>D c 



1 2 

l A,B 



e-*Q, y - -Q 2 - S a/3 S Q/3 



D a D a (j) + D c 



I 2 

l A,B 







And the junction conditions read 

3 



yu 
A 4 



y=0 



v=va 



y=VB 



K 

y 



K~ I -CTcSn + T. 



Cv 



K 



Bv 



(2.12) 



(2.13) 
(2.14) 
(2.15) 
(2.16) 

(2.17) 
(2.18) 
(2.19) 



The notation [X]_ indicates that we evaluate the quantity X on both sides of the 
brane and take difference, [X]1l = X + — X— 



3. The gradient expansion method 

In this section we derive the effective equations of motion using the low energy ex- 



pansion method (gradient expansion method) first proposed by Kanno and Soda [13| 
to study our scenario. In this method, the full five dimensional equations of motion 
are solved, at different orders, in the bulk by performing a perturbation expansion 
in the metric. The parameter of expansion is defined so that the low energy regime 
corresponds to a regime in which the energy density (p) on the brane is smaller then 
the brane tension (cr), p << a. In this regime, the parameter of expansion can be 
expressed as 

I ' 
L 



(3.1) 
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where I is the bulk curvature scale of the AdSs and L is the brane curvature scale. 
According to the parameter of expansion ( |3.1|) , the quantities are expanded as 

E„" = E<?> + E«" + Ej?> + • • • . (3.2) 

Then, the iteration scheme consists in writing the metric as a sum of local tensors 
built out of the induced metric on the brane, 



9^(y, ^) = a 2 {y, x) [V(^) + 9$(V, ^) + 9$(v, O + • ■ •] 



gfJiy 



y,x' 



n 



1,2,3,... 



(3.3) 
(3.4) 



Here, y is a generic point in which the Dirichlet condition is taken on the brane. 

In the following, we derive the zeroth order and first order solutions using the 
above scheme. 



3.1 Zeroth order solutions 

At zeroth order matter is neglected, we intend vacuum brane, and going at higher 
orders means we are considering perturbation of the vacuum solution as matter is 
added to the brane. The equations to solve are 



° n (0)2 v (0)Qy,(0)/3 
-£4 ^f3 - p 





12 



A,B 



4 q(J) _ Iq(°) 2 _ S (0)a/3 E (0 j 



/ 2 

A,B 



The junction conditions are given by 



E (o), _ ^ g (0) 



S (0), _ ^ g(0 ) 

E (o> _ 3 s;Q{0) 



v=va 



V=Vb 



(3.5) 
(3.6) 

(3.7) 

(3.8) 

(3.9) 
(3.10) 
(3.11) 



Integrating equation ( p.5[ ) and using the constraint ( p78|) we obtain the solution of 
the traceless part of the extrinsic curvature at zeroth order 



y(0)u 







And the trace part of the extrinsic curvature 



Q 



(0) 



4 

r , Q 



(0) 



IB 



(3.12) 



(3.13) 
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where we have inserted ( 3.12 ) into (|3.6| ) and " |±" denotes the solutions of each region 
of the bulk space-time. Inserting (|3.12f) and (|3.13| ) into (|2.12|) , the extrinsic curvature 
at zeroth order is 



JO) 



, • (3-14) 

± La,B ± 

In this order, the equation ( 3.14 ) gives the evolution of the extrinsic curvature in 
different regions of the bulk space-time which correspond to two different values of 
the AdS curvature scales, I a and Ib- 

From definition of the extrinsic curvature and using equation (|3.14j) we need to 
solve two different solutions of the metric 



1 d 



(o) 



~2dy 9,w 
1 d 



2dy 



q (o) 



e' 

T A 

e' p 

Tl 



7 HuV 



JO) 



< y < va 



y A <y <Vb ■ 



Integrating equations Q3.15| ) and ( |3.16j ) we get the zeroth order metric as 



ds 2 



dy 2 + a 2 (y, x)hf u (x)dx fJ 'dx 



where 



a_(y,x) 
a + (y,x) 



~d A (y,x)/l A 
d B {y,x)/l B 



(3.15) 
(3.16) 

(3.17) 



(3.18) 
(3.19) 



Here, we have integrated from a generic point y' to y such that (1a.,b — f V , e^ y ' x ^dy 
is the proper distance between a generic point and y. The tensor h^ u is the in- 
duced metric tensor depending on the brane coordinates. In fact, the factor a± is a 
conformal factor that relates the metric on the branes. In this case we have 



h c ~ 

h C+ 



a 2 h B 



where /irvf 



n C— brane uA 
iJau 5 ' i ">' 



(3.20) 

■ !>,. ■•-' , in' ■ (3-21) 

., . . ,,,, g^u brane an< i = g^ brane are the induced metrics on 

the C-brane, A-brane and B-brane, respectively. 

From the junction conditions (|3.9| ) and ( |3.11| ) we have the fine tuning conditions 
for A-brane tension and B-brane tension respectively, 

6 

U' 

6 

~Ib ' 

Furthermore, the junction condition (|3.10p yields 

o 3 / a — 1 



n 2 a A 



n 2 a B 



(3.22) 
(3.23) 



k a c = - 



a 



(3.24) 



-7- 



where we have defined a = Ib/Ia- This equation is a fine tuning condition for C- 
brane tension in terms of the curvature scales I a and Is of the slices of AdS^ bulk. 
The equation ( |3.22| ) - ( |3.24| ) also implies the relation for the brane tensions: 



a A + 2<7c + a B = . 



(3.25) 



The RS I model is obtained for a = 1 where C-brane is absent, <Jc = 0. For a < 1 
we have o~c > which correspond to inflation C-brane [16]. Various brane world 
models can be recovered by using relations ( |3.24|) and (|3.25| ) [17[ |18 , |19| , [20|] . 



3.2 First order solutions 



The aim of this subsection is now to solve the four basic equations (|2.13|) - (|2.16| ) at 
first order. In this order, the solution can be obtained by taking into account the 
terms neglected at the zeroth order. We have 



6 ^ U> 



-Q« = [^R] 



+±6;(D«D a 4> + D< 



(1) 



''AjB 



A,B 



\D a D a <p + D a (f)D c 



DJ§» ~ \D,QV = 



(3.26) 
(3.27) 
(3.28) 
(3.29) 



where the superscript (1) represents the order of the gradient expansion. The junction 
conditions are given by 



SW" _ ^g(i) 



— T 



y=0 



y=VB 



Au 

fj, > 



K T. 



K rpBis 
2 " 



(3.30) 
(3.31) 
(3.32) 



Here, T Au , T^ u and T^ u are the energy- momentum tensors with the indices raised by 
the induced metric on the A-brane, C-brane and B-brane, respectively. We compute 
the Ricci tensor [( 4 ).R^] ^ of a metric a^h^. The Christoffel symbol is 



KM = KM) + D^aSt + £Un< - h^D" In a . 



(3.33) 
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Here is a covariant derivative with respect to the metric h^. Using equation 
( |3.33|) the Ricci tensor of a metric g^ v = a 2 _h~ u is given by 



[ {4) R;(g)] {1) = \ \ {4) R;(h-) + ^(v,V»d A + ±V,d A V»d A ) 



+-8 u JV a V a d A - -V a d A V a d A 



Ia 



(3.34) 



Contracting indices fi and v of the equation ( ^.34[ ), we obtain the expression for the 
Ricci scalar 



(i) 



(4), 



6 



1 



R{h~) + t V oV' T d A - —V a d A V a d 



On the other hand, for the Ricci tensor and Ricci scalar of a metric 
respectively are given by 

+ [ ifi \ <*B 



(3.35) 
a 2 h+ 



I 



B 



--5»[V a V°d B + -V a d B V°d 



B 



[ (4) *rf = 4- 



6 



B 



R(h + ) - — V a V°d B + —V a d B ?rd B 



B 



(3.36) 
(3.37) 



We also express the kinetic terms of in terms of the proper distance, the indices 
— and + have been omitted to give a general CclSC, clS follows 

1 / 

v^v v <t> - v^v v <t> + -x ( v a v a <p + v a <pv a 



d 



a? dy 



4 m 



~| P tl <lV v d-~6r.V„dV a (l 



I 



4 c " 



(.3.38) 



In the equations above, Z> M denotes the covariant derivative with respect to the 
induced metric h^ u on the brane. And we have written the Ricci tensor in terms of 
the proper distance for a generic point on the y-axis. 

Substituting (|3.35|) into (|3.27|) we obtain the solutions of negative side of the 
C-brane as 



Q 



Wi 



Ia 

„2 



V a V°d A - —V a d A V°d A 



(3.39) 



We are now to obtain the traceless part of the extrinsic curvature at first order Y^ v . 
Substituting ( |3.35|) and ( [3.381 ) into (|3.26|) and integrating we obtain 

Ia 



1 



l A at 



V^d A V v d A - -Sp a dAD a d A 



^4 ' 



(3.40) 
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where x^\- is a integration constant which satisfy X^l- = an d — 0- 

And the solutions of positive side of the C-brane are 



Ib_ 

2 



a 



6 



B 



- (4) i2 - — ( V a V a d B - —V a d B V a d 



D 



B 



(3.41) 



B 



2a\ 
1 

IbcA- 



" 4 / 
1 

T 



V^d B V V d B - -S'PvdBWdB ) + 



(3.42) 



"B 



0. 



where is a integration constant which satisfy Xu\+ = ® an d 

The integration constants an d are non-local terms , which corresponds 
to the projection on the brane of the five- dimensional Weyl tensor. Therefore they 
carry the information of the bulk gravitational fields. In the following subsection we 
give explicitly the expression of these terms. 



4. The effective equations of motion on the branes 

In this section, we derive the effective equations of motion on each brane, at the 
first order, by substituting the equations of the extrinsic curvature into the junction 
conditions ( |3.31|) -( j02"|) . After imposing the trace of the projective Weyl tensor 
vanishes, Xu = 0, we then obtain the equations of motion for the scalar (radion) 
fields. We start to derive the equations of motion on the middle brane. 

4.1 Solutions on the middle brane 

At the middle brane (C-brane) is no Z 2 symmetry, the traceless part of the extrinsic 
curvature is asymmetric. Therefore it has two values in different sides of C-brane. 
Using equation ( |3.39| ) - (|3.42p , the junction condition at C-brane is written as 



, , , , — (a - 1) G u Jh c ) = k 2 T^ u . (4.1) 

e 4d B /l A e -4d A /l A 2 M M 

where we have used the Einstein equations on the negative and positive sides at 
C-brane, G u \- = G"\+ = G^(h c ). The junction condition at A-brane is given by 



+ ( V ^d A + \5;V a d A V°d^ = ^fcTf ,(4.2) 

and the junction condition at B-brane yields 

x;\ + + ^e 2d ^G^h c ) + e 2dB ' l ° (Vjrd B - 5;V a V°d B ) 

+j-/ dBllB (r>»d B v v d B + ^v^bV^b) = -!L e MB/iB T Bv ; (43) 
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where T>^ is the derivative covariant with respect to the induced metric on the C- 
brane. To get equations ( f4.2|) and (|4.3j ), we have used the conformal transformation 
of the metric both at negative side of C-brane and A-brane and at positive side of 
C-brane and B-brane as follows 



h c ~ 

h c+ 



-2d A /l Ah A 

c U flU 



(4.4) 
(4.5) 



Therefore, the index of T^ u and T® v are the energy-momentum tensors with the 
indices raised by the induced metric on the both side of C-brane, while T^ u and T^ v 
are he energy-momentum tensors with the indices raised by the induced metric on 
the A-brane and B-brane, respectively. In order to obtain the effective equations 
of motion at C-brane, one can subtract equation ( [4.2|) with respect to (|4.3|) . Then, 



substituting this result into equation ( fOD yields 

2k 2 



G»(h C ) 



l A {$c + W c ) 
1 

+ ($ c + a^ c ) 
, ^($c 



T C / + ±(1 + *c)Tf + \{l - *c)T* v 



c 



a 



c 



where we have defined two scalar fields 



(4.6) 



c 



2d,A/lA 



to>(<3> 



c 



c 



^> c = l-e 
\ 3 ^>c 



2d B /l B 



C 



Inserting 



Xu, I — 



X l 



2(l + $c) v ' 2(l + ^ c 

into ( |4.2| ) and ((4.3|), respectively, we obtain 
k 2 (1 + $ 



f c 









Ia 






l-*c) 


(<*>c 






al A 



T° v + i(l - a* )I)f + ^(1 - *c)T,f 



T- + i(l + ^)T- + i(^ 



T 



•c 



(4.7) 
(4.8) 



•c)] , (4-9) 



(4.10) 
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where 



c) 



c 



c 



c 



, V U V C V U V C 
ty c \ 2 



'c 



(4.11) 



(4.12) 



The equations ( |4.9|) and ( [4.1 0|) correspond with discontinuity of the evolution for 
Weyl tensor in each regions. The value Xul- on the negative side corresponds to the 
evolution toward A-brane and the value Xu\+ on the positive side corresponds to the 
evolution toward B-brane. 

The effective equations of motion for both scalar fields $c and can be ob- 



tained by using the condition Xu\- = anc ^ Xt 



0, respectively. Then we get 



c 



c 



K 



U(l-a) 
1 



1 - a)T A + 2T C 
2cu($ c ) + 3 
(iu;($c N 



2w($ c ) 

.2 



ft 



a/^a — 1, 
1 



3 d$c 
'{a- 1)T B + 2aT c 



c 



c 



3 dm 



c 



(4.13) 



(4.14) 



We see that the effective equations of motion for the scalar fields dependent on the 
matter sources for two branes. 



4.2 Solutions on the other branes 

We now proceed to find the equations of motion on the other branes. First, we derive 
the equations of motion on the A-brane. Using the fact that the metric induced on 
the A-brane is related by hA v = e 2dAlA h^~ , the junction condition on the A-brane is 
given by 



The junction condition at C-brane is 



(4.15) 



Ad B /h 



-Ad A /l a 
2d A /l A 



+ (a- l)e 



K 2 e 2d A /l A rpCv 



(V M V^ A - S?yS7 ff d A ) + - V^d A V v d A + -5y a d A V a d A 



(4.16) 
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and the junction condition at B-brane yields 

°^A 2d.A/lA+2d a /l. ™ /iA\ 



xl\+ + 



I e ^d A /l A +2d B /I b qv^ 



+ae 



2d A /l A +2d B /l B 



+e 



2d A /l A +2d B /l 1 



+ _ e 2d A /i A+ 2d B /i B ( y dA v dB + ^VvdAVds 



Ia 



1 

2 -/i 



— —. p1d A /l A +2d B /l BT Bv 

2 M ' 



(4.17) 



where V is a derivative covariant with respect to the induced metric on the A-brane. 
Eliminating X^l- an -d t ne equations of motion on the A-brane is given by 



2k 2 



Z A ($ A + a(l-$ A )* A ) 



l T ^ + ( l_^ ) ( T ov + ^ (1 _ lIrA)7 , 



2 ($ A + a(l - $a)®a) 
+aP;(<S> A ,V A ) 



a(l - *a) Wa) + (1 - «*a) Wa) 



(4.18) 



where 



^(*a) 



+ 2(i-$ A ) ( V ^ V ^ - ^; v ^^ vff ^ 

3 

+ 2(1-* A ) 



p;($ a ,^ a ) = v M $ A v* A + ~6»v a $ A w A . 

Here, we have defined two scalar fields as follows 

$ A = 1 - e ~ 2dA/lA , * A = 1 - e~ 2dB/lB . 
Substituting ( |4.18|) into (|4.15|) and (|4.16|) , respectively, we find 



(4.19) 

(4.20) 
(4.21) 

(4.22) 



X/i I 



K 2 (l-$ A )(l-gy A ) 

' ($ A + a(l - $ A )tf A ) 
Ia 



1 



_rp A u I 

2 " ^(l-a9 A ) 
a(l-$ A )P^ A ) 



t^ + ^(i-* a )t; 



i 



Bv 



2($ A + a(l-$ A )^ A ) 

+ (i - «^ A ) p;($ a ) + «p;($a, 



(4.23) 



13 



an 



Z ($a + «(1 - $ a)^a) (1 - - $ a) 

1 a(l - $4) + $a t bJ 



l T f + (1 - $4)^ 



cv 



+- 



a 



Oil A 



2 ($ A + a(l - $> a )Va) (1 - - $a) 



1 



(1-*a) 



(4.24) 



Using the conditions x%\- — and — we find the effective equations of 
motions for the scalar fields $a and ^a as follows 



2k 2 (1-$ a ) 



-{l-a)T A + T c 



Ia(1-<x) 



1 /a- 1 



a 



iB 



(4.25) 
(.4.26) 



Second, we now derive the equations of motion on the B-brane. Using the relation 
of the metric between C-brane and B-brane the junction conditions at the A-brane 
is 

Ia 



x ^_ + >A-2 dA /l A -2 dB /l BG ^ hB) 



_l_A -2d A /l A -2d B /l B 

Ib 



v„vdB - s;w a w a d B ) - — v^d B v u d B + -5;v a d B v a d B 



j^_ e -2d A /l A -2d B /l B 



1 

T A 



V^ a V^a + -S;V a d A V a d A 



+l e - 2d ^ lA - 2dB / lB (v»d A V u d B + Uft a d A V°d B ^ 

K 2 

- _„-2d A /l A -2d B /l BT Au 

2 " ' 

and the junction condition at the C-brane is given by 



(4.27) 



e id B /l A e -4d A /l 

2e~ 2dB ' lB [Ib_Ia 
= ^ e -^B/i BT Cu ^ 



_ I p -2d B /l B fl_B _ l_A_\ r u (h B\ 

-4d A /u + 6 y 2 2 ) ^ ' 



v M Vd B - sy a v a d B ) - — v^d B v u d B + -5 v l y tr d B v a d B 



(4.28) 



where V is the covariant derivative with respect to B-brane, while the junction 
condition at B-brane is 



al 



(4.29) 
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Eliminating Xu\- anc ^ from equation fl4.27Q -( |l~2S ), the equations of motion 
on the B-brane is given by 



Gl(h B ) 



2k 7 



where 



+ -(^ B + 1)($ B + 1)T / 



Ay 
ft 



(<f> B (q B + 1) + a* s ) 



(^ + i)p;($ B ) + ($ B + a)p;(^) 



3 



2($ B + 1) V ' "/x v --i,x 



P»{$ B , *b) = V^bV^b + -5;V ff $BV ff *B 



(4.30) 



(4.31) 

(4.32) 
(4.33) 



The two scalar fields are defined as follows 



B 



,2d A /l A 



^ B = e 2ds/ ' s - 1 . 



(4.34) 



The solutions for x^\- and Xa\+ can be obtained by substituting equation ( |4.30| ) into 
( |4.27| ) and ( |4.29| ), respectively. Then, we get 



K 



($ B (^B + 1) + a^ B ){^ B + 1)(* B + 1) 



/4 



P, 



At/ 



2 ($ B (^B + 1) + Qt* fl ) ($ B + 1) 2 (^B + I) 2 L 

+ ((* fl + 1) - a* B ) ((* fl + i)p;($ fl ) - p;($b, 

/t 2 (^ B + l)(q + $ B ) 
($b(^b + 1) + aV B ) 



a($ B + 1)P:(*b) 



(4.35) 



^-rp B u _|_ 


/ „ N 


2 " 


Va + $B/ 



p 



2 la + $ B 



A;/ 
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oJa 

"2($ B (^ B + 1) + aV B ) 



+ 1)P;($ S ) + ($B + B ) 



(4.36) 



Finally, the equations of motion for the scalar fields are given by 



1 



2k 2 ($ 



1) 



2/t 2 (^ + l) 
U(a-l) 



Ia 

rpB _|_ 



2 (1-a) 



(4.37) 
(4.38) 



In the derivation of equations of motion above we first to know the dynamics on 
one brane. Then we know the gravity on the other branes. Since the dynamics on 
each branes are not independent, the transformation rules for the scalar fields are 
given by 



c 



$R = 



1 - $, 



c 



ty A = 



(4.39) 



In the following subsection, for the realization at the first order expansion, we 
study the cosmological consequence of the radion dynamics. We need to derive the 
Friedmann equation on each branes. 

4.3 The effective Friedmann equation on the branes 

The metric induced on the brane is the Friedmann-Robertson- Walker (FRW) metric, 



ds 2 = -dt 2 + a 2 (t)'j mn dx m dx n , 



A,B,C 



where the time and space components of the Einstein equations are given by 

'of + k 



GT = 



laid; + a? + k 



(4.40) 

(4.41) 
(4.42) 



and k is the spatial curvature, k = 0, ±1. 

Assuming that the matter of the energy- momentum tensor to be = —p l 5^, 
(i = A, B, C), then the field equations on the middle brane can be written as 



G° (h c ) 



2JC 



2«y 



l A (<S>c + a^>c) 
3 



2a A 



2a B 



+ 



(<S>c + a^c) 



4\(l + $ c ) (l-* c ) 



(4.43) 
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2JJ 



fir 



2k 2 P ' 



1 

3{ & c 



;i + $ c ) (i-* ) 



20:^ 2ctB 

2H(<5> C + a^c) + ($c + a^ c 



1 + $ 



c, 



where the equations of motion for the radion fields are 



n 2 p c {l + $ 



ty c + 3#^ 



c 



3U 

K 2 p C {l - ^ 



1 



3L 



1 



2« A 



+ 



2(1 + $ 



c 



[a 



2a 



B 



2(1 - * 



(4.44) 

(4.45) 
(4.46) 

(4.47) 

Using the equations (|4.41 ) and ( 4.42 ) and eliminating <E> and ^ from ( 4.45|) and 
( [4.46| ), respectively, we get 

k An 2 



c 



Here, we have defined two dimensionaless parameters 



a A 



P 



c 



OLE 



P 



C 



,B 



H + 2H + 



3l A (a — 1 



-P 



c 



(4.48) 



Integrating this equation, we obtain the Friedmann equation with dark radiation 

k 2k 2 

—n 

a 



H 2 + 



C c 



(4.49) 



a 2 -; 3l A (a — 1) 

where Cc is an integration constant. 

We derive the Friedmann equation on the A-brane where the FRW metric is given 
by ( |4.40| ). From the equation ( |4.18| ) the time component of the Einstein equation is 



OruA\ 



G»(h 



+ 



K 2 p C 



w 

Ia 



2a A 



2a 



B 



219 



a(l - <t> A )P°(V A ) + (1 - a* a) P$(®a) + *a) (4-50) 



where 19 is defined as follows 



W 



($ A + a(l-$ A )V A ) 



(4.51) 



And the time component of the kinetic terms are given by 

3 



Pg(* A ) = 3H* A 
P °($4) = 3H$> A 



4(1 ^ A, 

3 

4(1-$,) 



^ A j 


(4.52) 


* A ' 


(4.53) 




(4.54) 
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The space component of the Einstein equation ( 4.18 ) is given by 



G rn {h A, 



K 2 p C 
W 

Ia ' 
216 L 



2a A 



(1 - $ 



A 



2a B 



a(l - $ + (1 - a* A ) iC(^) + *a) (4-55) 



where 



P^a) 
P^a) 
P^a^a) 



+ 2H^ A + 



$ A + 2H$ A + 



4(1 -* A ) 
3 

4(1 - $a) 



* A i 



^ A i 



(4.56) 
(4.57) 
(4.58) 



Inserting equations ( |4.52|) -( [4.54|) into (|4.5Q) and (|4.56| )- p34]) into ( f4.55|) , respec- 
tively, we obtain 



3 [H 2 + 



k 



''A 



K 2 p C 



+ 



w 

Ia 



+ (1 + 77^(1-^) 



2« A 

a(l-$ A )^3^ A - 



2a s 
3 



2Z(9 

+(1 - a* A ) ( 3#$a 



4(1 
3 



.gr- 



and 



2LH"-t -3 ^+-?- 



k 



K 2 p C 



4(1 - $ A 



3 • 

2$ A 



(4.59) 



16 



2a, 



+ {1-$ A ) 1 + -— 



2a 



+ 



2Z(9 



a(l - 



a; i ^a 



4(1 -^a 



-*a 



-$a*a 



(4.60) 



4(1-$ A )~ A , 

where the equations of motion for the radion fields $a and ty A are obtained from 
the equations (|4.25| ) and (|4.26|) , respectively 



$a 
*A 



%K 2 p c (l-§ A ) \{l-a) 
2>1 A {1 — a) |_ 2oa 

8k 2 p C (1 - 
3/ A (l-«) 

+ — 1 , > A^A • 



1 - 



(i-q) 

2a«s 



3i/g, - 



2(1 - $4 
1 



2(1 -*a) 



^ A i 



(4.61) 



2(1 -$a) 

Substituting ( |4.61| ) and ( |4.62j ) into (|4.60|) , respectively, we obtain 



2„A 



H + 2H Z + ^- 



k 2n z p 



3/a 



(4.62) 



(4.63) 
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Then, we obtain the Friedmann equation with dark radiation by integrating ( 4.77 ) 
on the A-brane as follows 



H 2 + ^ 



«V | C A 

31a cl a 



4 ' 
A 



(4.64) 



where Ca is an integration constant. 

Finally, To obtain the Friedmann equation on the B-brane we use the same 
procedure. The Einstein equations on B-brane are given as follows 



-3 ( H 2 + A 

a 2 B 



K 2 P C 



+ 



w 

Ia 



2a B 

1 + ^ B )[3H^ B + 



+ (1 + *b) 1 + — (1 + $ 



2a A 



219 



4(1 + $ B 
3 



-cb 2 



4(1 + * B ) 



3 ■ 



*B + 



(4.65) 



and 



2 # 



3 H 2 + 



"B 



K 2 p C 
10 



2a 



B 



+ (1 + *b) 1 + ^ — + * 



2«A 



+ - 



1 + *b) $B + 2if$ 



219 



+ (a + <S> B )[ * B + 2F* 



B 



4(1 + $ s 
3 



-$1 



4(1 + 



*B 



2®aVa 



(4.66) 



where Z# is defined as follows 



Ia 



16 = —($ B (y B + l) + ay B ) 



(4.67) 



In order to obtain (|4.65 ) and ( |4.66| ) we have used the components of the scalar kinetic 
terms: 



P °(^ B ) = 3HV B + 



P °($ B ) = 3H$> B + 



P^b^b) 



4(1 + ^ B 
3 

4(1 + 



-\i/ 2 

* B ) 



B 5 



$B^B, 



p;(f B ) = f B + 2Ff B 



P" 1 ^ b) = $B + 2# <£> B - 



4(1 + ^ B 
3 

4(1 + $b) 



.\i/ 2 

* B ) 



(4.68) 

(4.69) 

(4.70) 
(4.71) 

(4.72) 

(4.73) 
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The equations of motion for the scalar fields $b and ^b are used to eliminate 
the second derivative <3>s and in the equation ( |4.66[) , where it is given by 

8k 2 P c (1 + * b ) 



31a(cx — 1 
Ik 2 P C (1 + $ 



B 



3l A 

1 



2aa.B 
1 



lOLA 



+ 



+ 1 
1 



a) 



3H^ B + 
-3H&B 



2(1 + 

1 



(4.74) 



2(1 + $ 



B 



2(1 + » B )-'"- (475) 

Solving the equation ( |4.6(j| ) we get the Friedmann equation on the B-brane by inte- 
grating the equation below 



H + 2H 2 + X- 



to find 



H 2 + 



n 2 p B 



2k 2 P b 
3l B 

C b 



31 



+ 



B 



J 4 ' 
1 B 



(4.76) 
(4.77) 



where Cb is an integration constant. 



5. The scalar-tensor gravity 

In the previous section we have derived the effective equations of motion in a three 
brane system. Now we show how we can write the scalar-tensor gravity using the 
effective equations of motion on this system. We use the solutions on the middle 
brane to obtain a scalar-tensor gravity with two independent scalar fields. In the 
following we omit subscript C of the equations that related to the middle brane. 
From the equation ( |4.6| ) we see that a term ^ ($ + a^/) can be defined as a first 
dimensionless scalar field, 

/0^k($ + a vD) , (5.1) 

where I is an arbitrary unit of length. Because the scalar fields $ and ^ correspond 
to the proper distance, the definition of the scalar field (|5.1[ ) associated with overall 
distance of the middle brane. Then, the second scalar field can also defined as a 
function of both scalar fields, we define 

= . (5.2) 



We intend to write the effective equations of motion on the middle brane (|4T 
as follows 

= /! + $ Av 1-tt 
" l(j) V 2 " 2 



(5.3) 
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where u(4>) and u(<f>) are the arbitrary functional coupling of 0. The absence of 
mixing terms in the equation above yields the following constraints 



21, 

Applying this constraint into equation 
- 2 x l + $„„,. 1 - * 



5$ «^ 
we have 







rpAv I 

2 



7? " + if) + i (2?^ - ^D0) 



^(0) 



2/ 



-u(<j>) 



uJ(0) 



On the other hand, by inserting ( |5.1| ) into ( [4.6| ) we obtain 



K 2 ( 1 + $ 



/0 



rpAu 

2 M 



rpBv _i_ rpCv 

2 M ' ^ 



(P^0 - 



3/ 



.4 



4/0(1 + $) 



4/0(1 - 

The relation between the coefficients in equation ( |5.5|) and ( |5.6|) is 



3l A 



4/0(1 + $) 



^(0) (l_A 
21 



21 



a 



d£/d<$> 



and using the constraint (^) we get an equation of the form 



;i + $) d£/d<$> 



a — 



d£/d<$> 



{l-ty)d£/d^ 

It is easy to see that the solution of the equation ( |5.8j ) are 



. 



ft- 



(1 + $) 

i - m . 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



(5.8) 

(5.9) 
(5.10) 



The first solution yields vanishing the coefficient of ( |5.5|) and ( |5.6|) . Then we find a 
solution, 



(5.11) 
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Substituting this solution into (|5.4j) we obtain a differential equation for p 



dp _ y^eg/ 2 {{l-a)l A /2 + l(j>f 
di ~ (e« - a) V Z 2 ^ 

In order for to be only a function of £, we require that in equation ( p,12[ ) 

' (1 -^+^V=^)- 



v(0 



(5.12) 



(5.13) 



Then, we obtain the solution of tp, 



= a coth 2 



(5.14) 



Finally, the effective equations of motion on the middle brane can be written as 



^"/0 



3/0 



my 



cosh 2 (^)T^ + 



3/0 



_2l A u(cj)) V 2 J M 2/ b oj(0)~ V 2 



sinlrf ^ l^+Tf" 



4a; (0) 



where 



^(0) 



Z0 



(5.15) 



(5.16) 



2 V^+(l-a)W 2 , 
The effective action for the the middle brane corresponding to the effective equa- 
tions of motion (15. II) can be rewrite as 



S = — - / d 4 xV—h 
2k 2 



uU) 9 



+ / d A xV^h 



£ c + -J^-cosh 2 ^ 



3/0 



2 fVN /.B 



, , _ , sinh — ) £ 

2/^(0) V27 2/ b cj(0) V2/ 



(5-17) 



where C c and £ B are the Lagrangian correspond to A-brane, C-brane and B- 
brane, respectively. This action is the scalar-tensor gravity on the brane with two 
scalar fields as a function of the two proper distance. 



6. Conclusions 

In this paper, we consider three 3-brane systems with with A and B-brane are placed 
at the fixed point of the orbifold whereas the C-brane is put between A- and B-brane. 

We use the gradient expansion method to analyze, in the first order, the effective 
equations of motion, in particular the radion Lagrangian. In this case we derived the 
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Friedmann equation with dark energy radiation by direct elimination of the radion 
fields in the Einstein equations. We also derive the scalar-tensor gravity with depend 
on the scalar (radion) fields. 

We can also generalize this scenario to the multi (more than three) 3-branes 
system in the low energy limit. It is also interesting, in this model, to investigate for 
higher order correction. 
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